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Abstract. We study the secular instability of magnetized differentially rotating radiative zones taking account of viscosity and 
magnetic and thermal diffusivities. The considered instability generalizes the well-known Goldreich-Schubert-Fricke instability 
for the case of a sufficiently strong nagnetic field. In magnetized stars, instability can lead to a formation of non-spherical 
unstable zones where weak turbulence mixes the material between the surface and interiors. Such unstable zones can manifest 
themselves by a non-spherical distribution of abundance anormalies on the stellar surface. 
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1. Introduction 

The effect of the magnetic field on stellar differential rota- 
tion was considered by many authors. In the pioneering work, 
Fricke (1969) first pointed out the importance of thermal dif- 
fusion for stability of magnetized differentially rotating stars. 
This conclusion had a futher development in the paper by 
Acheson (1978) who took into account consistently the main 
diffusive processes (viscosity, thermal and magnetic diffusion) 
and noted that even such a weak kinetic effect as viscosity 
can influence the secular stability of magnetized stellar radia- 
tive zones. Acheson (1978) reviewed in detail the so called 
double diffusive instabilities that can occur in stellar radiative 
zones because of dissipative effects. Perhaps the best known 
in astrophysical context diffusive instability is that considered 
by Goldreich & Schubert (1967) and Fricke (1968). In non- 
magnetic stars, this instability arises if the angular velocity de- 
pends on the height above the equatorial plane and can be re- 
sponsible, for example, for the mixing and angular momentum 
transport in stellar radiative zones (see, e.g., Knobloch & Spruit 
1982, Korycansky 1991). Note that some double diffusive in- 
stabilities can be of importance for the rotational evolution of 
degenerate stars as well (Urpin 2003a, 2003b). 

Most likely, however, that the magnetic field may exist in 
stellar radiative zones despite many uncertainties regarding the 
origin, evolution and geometry of these fields. It seems that 
even magnetic buoyancy cannot substantially reduce the field 
in radiative zones because of a stabilizing influence of stellar 
rotation on magnetic buoyancy instability. Generally, rotation 
can suppress this instability if the rotation velocity, V, is larger 
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than the Alfven velocity, V > ca (see, e.g., Gilman 1970, 
Acheson & Gibson 1978, Acheson 1979). 

Secular instabilities are typically sensitive to the magnetic 
field, and the stability properties of magnetic and non-magnetic 
secular modes can be substantially different. The influence of 
the toroidal magnetic field on the Goldreich-Schubert-Fricke 
instability has been considered by Acheson (1978) who noted 
that such a field plays a destabilizing role and can change qual- 
itatively the criterion of rotational instability. This effect can 
be of importance, for example, for the sun where the toroidal 
field is expected to be buried deep inside the radiative zone 
(Schmitt & Rosner 1983). Note that a destabilizing influence of 
the toroidal field is not surprising because such field leads to in- 
stability even in non-rotating stars (see Tayler 1973, Goossens 
& Tayler 1980). 

The effect of the poloidal field can generally be more com- 
plicated and, to the best of our knowledge, has not been con- 
sidered in its full generality. Recently, Menou, Balbus & Spruit 
(2004) considered diffusive instabilities that occur in magne- 
tized differentially rotating stellar radiative zones and obtained 
the criteria of instability in the limit when either viscosity or 
magnetic diffusivity is vanishing. The consideration, however, 
was restricted by the case of a very weak magnetic field. The 
authors have found that stability properties are sensitive to the 
combination of viscosity, magnetic and thermal diffusivities 
and, following Acheson (1978), concluded that all three should 
generally be included in the analysis from the very beginning. 
The results of this consideration have been applied by Menou 
(2004) to account for millisecond oscillations in differentially 
rotating atmospheres of accreting weakly-magnetized neutron 
stars. 



Ihe present paper extends the study ot secular instability to 
the case of a relatively strong magnetic field and arbitrary dif- 
fusivities. We show that such the magnetic field changes dras- 
tically the stability properties of radiative zones. In a strong 
field, the instability of differential rotation can lead to a forma- 
tion of non-spherical unstable zones where weak turbulence is 
generated. The geometry of these zones depends on the mag- 
netic configuration and rotation. For example, in the simplest 
case of an axisymmetric field the turbulent zones are formed 
near the rotation axis and around the equator. If the field is not 
axisymmetric then the instability will lead to the formation of 
cone-like unstable zones within the radiative region. Turbulent 
motions mix the material in the unstable zones, thus a fraction 
of the material from deep interiors may reach the surface where 
chemically pecular regions should be formed. 

An outline of the paper is as follows. In Section, we de- 
rive the general dispersion equation taking account of viscos- 
ity, magnetic and thermal diffusivities. The criteria of instabil- 
ity are discussed in Section 3, and the growth time of secular 
instability is considered in Section 4. A discussion of results 
follows in Section 5. 

2. Basic equations and the dispersion relation 

The MHD equations governing the velocity, magnetic field and 
thermal balance in the stellar radiative zone are 



v + (v ■ V)v 



P 



1 

Airp 



(VxB)xB 



Av + ^V(V • v) 



p + V • (pv) = 0, 
B-Vx(bxB) 
V • B = 0, 



-Vx [ij.(VxB)], 



pc p f -f3p + pc p v ■ (AWT) = V • (k ■ VT) 



(1) 

(2) 
(3) 
(4) 
(5) 



(see, e.g., Landau & Lifshitz 1981). Our notation is as fol- 
lows: v is the fluid velocity; p, p, T are the density, pres- 
sure and temperature, respectively; B is the magnetic field, 
g is gravity, AVT = VT - V ad T is a difference of the 
real and adiabatic temperature gradients, V a dT = f3\7p/pc p ; 
(3 = — (d In p/d In T) p , c p is the specific heat at constant pres- 
sure; v is the kinematic viscosity; rj and k are the magnetic 
diffusivity and thermal conductivity, respectively. We ignore 
the spatial dependence of the kinetic coefficients v, rj, and k, 
which is appropriate for a short wavelength analysis. 

We work in cylindrical coordinates (s. ip, z) with the unit 
vectors (e s , e v , e z ) and assume that the fluid rotates with the 
angular velocity il — £l(s, z) in the unperturbed state. We as- 
sume the Alfven speed, ca = B/y/Anp, to be small compared 
to the sound speed, c s . Then, the magnetic field should satisfy 
the condition 



B<3x 10 7 a 



G, 



(6) 



where p is measured in g/cc and Tq = T/1Q 6 K. This assump- 
tion still leaves a wide range of B which provide an appreciable 



influence on the stability properties and are ot interest tor stel- 
lar radiative zones. Note that under the condition c s > ca, the 
Alfven speed can be larger than the rotation velocity sCl. 

In the unperturbed state, the star is assumed to be in hydro- 
static equilibrium in the s- and z-directions, 



— = G , G = g + Q 2 s + — (V x B) x B. 
p ' Anp 



(7) 



Generally, even if c s > ca and the unperturbed Lorentz force is 
small compared to the pressure force, the effect of the magnetic 
field cannot be neglected in this equation. Taking the curl of 
equation (7), we obtain 



■^Vj) x Vp 

P 2 

where 



on 2 r 



(8) 



L = 



A-k 



-V x 



-(V x B)xB 



Therefore, the magnetic field (as well as rotation) can be re- 
sponsible for the baroclinic effect that is oftenly crucial for sta- 
bility. The Lorentz force splits generally into potentional (mag- 
netic pressure) and solenoidal (magnetic tension) parts. In our 
study, the effect of the both parts on the unperturbed state is 
taken into account consistently. 

Consider stability of axisymmetric short wavelength per- 
turbations with the spacetime dependence oc exp(7i — ik ■ r). 
where k — (k s , 0, k z ) is the wavevector, \k-r \ S> 1. Small per- 
turbations will be indicated by subscript 1, whilst unperturbed 
quantities will have no subscript, except for indicating vector 
components. The linearized MHD-equations read in the local 
approximation 

(7 + w v )vi + ^k(k ■ ui) + 2fl x vi + se v (vi ■ Vfi) 

= ik— + G— — -r~— (fc x Bi) x B, (9) 

p p Airp 



7Pi + k ■ vi = 0, 

(7 + cu v )B 1 - se^Bx ■ Vfi) = -i(k ■ B)vi, 



k B x = 0, 



7/3 



(7 + w x )Ti - —pi = -vx ■ (AVT), 
pc p 



(10) 
(11) 



(12) 



(13) 



where the characteristic inverse timescales of the viscous and 
ohmic dissipation and the thermal diffusion are given by uj v = 
vk 2 , uj v — i]k 2 , and uj x — xk 2 , respectively; \ — K /pc p is the 
thermal diffusivity. 

Perturbations of the pressure in equation (9) can be ex- 
pressed in terms of perturbations of the density and tempera- 
ture, using the equation of state, 



pi = / dlnp \ /ft I\\ 
p \d\np) T \p^ P T J' 



(14) 



We consider the dispersion relation in the case kc s » 
7 that corresponds to the Boussinesq approximation for a 
non-magnetic fluid. Apart from this, the applicability of the 



Boussinesq approximation in the magnetic held imposes a 
restriction on the held strength. Namely, the magnetic held 
should satisfy the inequality c A /c 2 <C (kH)^ 1 where H is the 
lengthscale of unperturbed quantities. Solving equations (9)- 
(14) under these two conditions, we obtain the following dis- 
persion relation 

7 5 + a 4 7 4 + a 3 7 3 + a 2 -f 2 + ai7 + a = , (15) 
where 

<2 4 = 2bJ. q + 2uJ u + LO Xl 

a 3 = {uj, q + uj v ) 2 + 2{oj, n u v +LJ x u, n +Lj u Lj^ + 2Lj\+u 2 g + Q 2 , 
a 2 = oj x (co v + u n ) 2 + 2Wr)U) v (u) u + uj v + uj x )+ 

uj l/ (uj 2 + 2uj A ) + 2Lu v {uj 2 g +Lu 2 A + Q 2 )+Lu x (Q 2 + 2lo 2 a ), 
ai = 2(oj 1/ u) v + u) A ){u) u uu x + u u uu v + Ul) v uj x ) + 



u>n[{^-q + 2uj l y)uj 2 + {ur, + 2u> x )Q 2 



u\ (cj 2 +lu 2 a + F 2 ) , 
a = u 2 (u x loI + u u uj 2 + u x Q 2 ) + 2uJ a uj v uj 7} u) x + 
u\ (u n u 2 + w x u\ +uj x F 2 ) , 

and 
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k B 



k 



G--(k-G) 



C = — AVT . 
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This equation describes hve low-frequency modes that exist in 
differentially rotating and weakly magnetized non-ideal fluids. 
In the limit of ideal MHD (y = r\ = \ = 0)> equation (15) 
reduces to a quadric one that describes four dynamical modes. 
The hfth root corresponds to a secular mode that appears only 
due to dissipative effects, and we address the properties of this 
mode. In a vanishing magnetic held (B = 0, i] = 0), this 
mode has been considered by Goldreich & Schubert (1967) 
and Fricke (1968). Equation (15) differs only by notations from 
the dispersion equations derived by Menou, Balbus & Spruit 
(2004). Note, however, that these authors considered the sta- 
bility properties of Equation (15) only in a weak magnetic 
held, ft 3> u>a, and in either the inviscid or perfect-conductor 
limit. In this paper, the secular stability will be considered for 
a stronger magnetic held with uja > ^ (but still satisfying the 
condition c s > ca) and for arbitrary diffusivities 

3. Criteria of instability 

The Hurwitz theorem presents a sufficient condition of insta- 
bility and states that equation (15) has unstable solutions with 
a positive real part if one of the following inequalities 

a 4 < , a < , 
A\ = a 4 a 3 — a 2 < , 

A 2 = a 2 (a 4 a 3 - a 2 ) - a 4 (a 4 ai - a ) < , 

A 3 = (a 4 ai - a )[a2(a4 : a 3 - a 2 ) - a 4 (a 4 ai - a )]- 



-a (a 4 a 3 — a 2 ) 2 < 



is fulfilled (see Aleksandrov, Kolmogorov & Laurentiev 1985). 
Since w v , oj v , and oj x are positive dehned quantity, the condi- 
tion a 4 < never applies. 

We consider only the criterion ao < that corresponds 
to double diffusive instability and is a sufficient condition of 
this instability. We focus on relatively large wavelengths, A = 
2ir/k, for which oj a » w„w,. This condition can be rewritten 
as 



A > A c = 27r 



= 2tt- 



CA 



(17) 



where B p = B 2 + B 2 is the strength of the poloidal held. 
Mention that A should be small compared to H because we 
use a local approximation. Following Menou, Balbus & Spriut 
(2004), we can estimate put] ~ (2 — 6) x 10 3 g cm/s 2 in the ra- 
diative zone. Then, we have A c <~ 10 3 /B cm where B is mea- 
sured in Gauss. This leaves a comfortable range of wavelengths 
for which H > A > A c and our local analysis is valid. Note that 
although conditions for instability can indeed be found using 
our approach with the accuracy in terms <~ (A c /A) 2 , the critical 
wavelength that discriminates between stability and instability 
cannot be addressed by this analysis because this wavelength is 
<~ A c . 

Under the condition (17), the criterion ao < reads 



X (F 2 +uj 2 a )+ V lu 2 <0. 



(18) 



Equation (18) differs from the condition of instability obtained 
by Fricke (1968) by the last term on the r.h.s. Usually, rj/x is 
small in stars but, on the contrary, the ratio tu 2 /F 2 is large. 
Therefore, whether or not the contribution of the last term is 
small will depend crucially on the particular system under con- 
sideration. This point concerning viscosity and magnetic diffu- 
sivity has been hrst emphasized by Acheson (1978). The condi- 
tion (18) also differs essentially from the well-known criterion 
of the magnetorotational instability, 



F 2 + uj\ + uj 2 < 



(19) 



(e.g., Fricke 1969, Balbus & Hawley 1991, Balbus 1995, Urpin 
1996). Note that criterion (19) is valid only if dissipative ef- 
fects are negligible and the magnetic held is very weak, oja <C 
s|fi s ,z|- Compared to the magnetorotational instability, the sta- 
bilizing effect of stratification is much reduced in equation 
(18). Since usually cu 2 > F 2 in stars, the criterion (18) re- 
quires a much weaker differential rotation for instability than 
the criterion (19). 

In a realistic situation, all three terms in equation (18) 
should be taken into account. Generally, the buoyancy fre- 
quency is much larger than il, and the stratihcation term is 
not negligible even in a fully ionized stellar plasma where r] 
is small. The inequality (18) depends on the direction of a 
wavevector and can be rewritten as follows 

N (k)^F 2 + . 2 A + ^ 2 g = ^ 2 + A f 2 - 



(16) 



D- 



. k s k 2 



k 2 



k 2 

4<o. 



(20) 



where 

A = sfl 2 s + uj 2 AQ cos 2 a + -C z G Zl 
X 

D = sQ 2 - uj 2 A0 sin 2a - ^{C S G Z + C Z G S ) 7 
X 

E = u A0 sin 2 a + -C S G S . 

X 

In these expressions, we use notations 



C G , lo ao 



k 2 B 2 p 
4wp 



for the Brunt-Vaisala and Alfven frequences, respectively; a is 
the angle between the poloidal magnetic field and the axis of 
rotation. In what follows, we use the procedure proposed by 
Miralles, Pons & Urpin (2004) to derive the criterion that does 
not depend on the direction of k and determines the minimal 
differential rotation leading to instability. Since the dependence 
of N(k) on the direction of k is rather simple, we can obtain 
that minimum of N(k) corresponds to 



fc| = 1 
k 2 ~ 2 



1 ± 



{A-Ef 



(A - E) 2 + D 2 



(21) 



Then, the minimum value of N(k) yields the following condi- 
tion of instability 



sn 2 s + + u 2 A0 ± ^D 2 + (A - E) 2 < 0. (22) 

X 

The gradients of p, p and T are related by 
Vp _ fdlnp\ (V£ VT\ 



P 



dlnp 



+ 0- 



T J 



(23) 



T \ P 

It can be readily obtained from equation (8) that the condition 
of hydrostatic equilibrium leads to 

sn 2 z = [Cx G} v = C Z G S - C S G Z - V (24) 
Then, equation (22) can be further simplified to obtain 

sfll + + u 2 AQ ± | ^ 2 + ^ 2 + + 

4G z ^[s{C z n 2 -C s n 2 z ) + 0J 2 AQ {C z cos 2a+C s sin 2a)] 
X 

-4sin 2 <W^ ^-w 2 + s0 2 ^ -2^1--^ sJl 2 ^ sin 2a 



1- 



X 



sfl 2 — —Lu 

X 



1 2 



+ 2— L,„uo 2 An sin 2a 



1/2 



X 



< 0.(25) 



The two conditions for instability follow straightforward from 
the above expression: 

sfl 2 + ^-cj 2 + uj 2 A0 < 0; (26) 



X 



AG Z ^ [s{C z n 2 s -C s n 2 z ) + u 2 A0 (C z cos 2a+C s sin 2a)] 
X 

-4 sin 2 auj 2 A0 ( - uj 2 + sil 2 ] - 2 ( 1 - - ) sQ 2 z u 2 A0 sin 2a 



+ 



X 



+ 2-L^uj 2 m sin 2a > 0. (27) 

X 



Ihese two conditions look like the Solberg-H0iland condi- 
tions (Tassoul 2000), but with additional terms due to the mag- 
netic field and diffusive effects. Equations (26) and (27) are 
equivalent to the condition ao < with the accuracy in terms 
~ (A c /A) 2 . Therefore, the boundary between the stable and 
unstable regions can be determined by these conditions with 
the same accuracy. If B = and X] = 0, equation (27) re- 
duces to the well-known Goldreich-Schubert-Fricke criterion 
of instability, stating that differential rotation is always unsta- 
ble if VI 2 ^ (Goldreich & Schubert 1967, Fricke 1968). 
In the case of a non-rotating star (11 = 0) with B v » B p 
when the terms containing loao can be neglected, equation (27) 
yields rfL^/x 2 > that is always satisfied. Therefore, the 
star with a presumably toroidal field is always unstable that 
was first argued by Tayler (1973). Generally, if the poloidal 
and toroidal fields are more or less comparable we have L v <~ 
uJ A0 (kH)~ 2 <C uj ao , and the terms proportional L v can be 
neglected. 



3. 1. Instability in a weak magnetic field 

We consider equations (26) and (27) in the limiting cases of 
weak and relatively strong magnetic fields. We refer the mag- 
netic field weak if O 2 3> uj ao (it is assumed that sfi 2 ~ sQ 2 <~ 
fi 2 ). This condition can be written as 



Bp < B c 



6 x 10 2 VpA 8 fi_ 5 G 



(28) 



where A 8 = A/10 8 cm, 0_5 = £!/10~ 5 s, and pis measured in 
g/cc. The quantity B cr depends on p and, in general, the same 
field can be considered weak in the inner layers and strong in 
the outer layers. Note that the field should also satisfy the con- 
dition c s > ca (equation (6)) for validity of our analysis. It 
can be convenient to represent the inequality (28) in terms of 
the equatorial velocity, V — flR, where R is the stellar radius. 
Then, we have from equation (28) 



V > UBpaRnX^p-^ 2 km/s, 

where B p3 = B p /10 3 G, R n = i?/10 n cm. 
In a weak field, equations (26)-(27) yield 

sn 2 s + ^-uj 2 <0: 
X 

4GA s (c z n 2 s -c s n 2 z ) + ( s o 2 ) 2 > o. 

X \ X/ 



(29) 



(30) 
(31) 



The criteria turn out to be independent of the magnetic field in 
this limit and, under certain conditions, can be satisfied. 

Since > in the radiative zone, the inequality (30) 
can be satisfied only if £1 decreases with the cylindrical ra- 
dius. Taking into account that, typically, rj/x ~ 10~ 4 and 
uj 2 ~ 10~ 6 — 10~ 7 s~ 2 in the radiative zone and assuming 
that \sfl 2 \ <~ £1 2 , we obtain that the condition r]uj 2 /x < \sQ 2 s \ 
is satisfied if £1 > 10~ 5 w 3 s _1 where w 3 = ^o/10~ 3 s. 
Therefore, the instability caused by the condition (30) may oc- 
cur in relatively rapidly rotating stars with £1 > 10~ 5 wo3 s_1 
or V > 10i?iiw 3 km/s if dCt/ds < 0. 



Ihe condition (31) depends on the vertical gradient ol the 
angular velocity. If O 2 = then equation (31) yields 



4<4 C z G z Vl 2 s > 0. 
X 



(32) 



We have C Z G Z < since both C and G are approximately 
radial and oppositely directed in the radiative zone if the star is 
far from the rotational distortion. Therefore, the inequality (32) 
is equivalent to Q 2 < 0, and both slowly and rapidly rotating 
stars can be unstable in this case if the wavelength of perturba- 
tions satisfies the condition (17). 

If Q 2 . 7^ then the vertical dependence of the angular ve- 
locity can produce either stabilizing or destabilizing effect de- 
pending on Q(z). If rotation is slow and i]u; 2 /x > fi 2 , then the 
criterion (31) is equivalent to 



G z (c z n 2 -c s n 2 z ) >o. 



(33) 



In slowly rotating stars, G and C are approximately radial and 
directed inward and outward, respectivetely. Then, we have 



G z 



-Gcos0, C, «Csin<9, C z «Ccos( 



(34) 



where 9 is the polar angle. Substituting these expressions into 
equation (33), we obtain 



n 2 cos 2 9 - til sin0 cosf? < 0. 



(35) 



Hence, VL 2 provides a destabilizing influence if f2 increases 
with the distance from the equator and a stabilizing influence in 
the opposite case. If fi 2 < and f2 2 < then instability occurs 
only in a cone around the rotational axis, and the opening an- 
gle of this cone depends on the radial and vertical lengthscale 
of O. On the contrary, if fi 2 > and VL 2 > then the instability 
arises only in a region around the equatorial plane. 

Rapidly rotating stars with a sufficiently strong vertical dif- 
ferential rotation, \sQ 2 \ > t]uj 2 /x, can be unstable indepen- 
dently of the particular shape of Q(z) if A satisfies condition 
(17), and the instability occurs everywhere within the radiative 
zone. 

3.2. Instability in a strong magnetic field 

We refer the magnetic field strong if 

B p > B cr = 6 x 10 2 VpA 8 f2_ 5 G (36) 

(but the condition (6) still holds) because such a field can pro- 
vide a strong influence on the stability properties. Equation (36) 
can be rewritten as 

V < 17 B^RuX^p- 1 / 2 km/s. (37) 
Under this condition, the criteria (26) and (27) transform into 

^ujI +uj 2 ao <0; (38) 
X 

2G Z -(C Z cos 2a + C s sin 2a) - 2 sin 2 a ( -u>l + sQ 2 s 
X V X 



- l 1 -^) s^ 2 sin 2a > 0.(39) 



Ihe first criterion never applies in the radiative zone since 
luq > 0, and we consider only the condition (39). Note that 
this criterion turns out to be independent of the strength of the 
magnetic field but depends on its direction. 

In slowly rotating stars with rjuo 2 /x > fi 2 , we have from 
equation (39) 



G Z (C Z cos 2a+C s sin 2a) — w 2 sin 2 a > 

or, using equations (34), 

(cos 9 cos a — sin 9 sin a) 2 < 



(40) 



(41) 



that never applies. Therefore, the radiative zones of slowly ro- 
tating stars with a strong magnetic field satisfying the inequal- 
ity (36) should be secularly stable. 

Consider the case of a relatively rapid rotation, £1 2 > 
rjoj^/x, or V > 10i?no;o3 km/s. Combining this condition 
with the inequality (37), we obtain the range of V where ro- 
tation can be considered rapid and the magnetic field strong 

17B p3 .RiiA^ V 1/2 km/s > V > 10i?nw 03 km/s. (42) 

Equation (39) yields within this range 

sin a(f2 2 sin a + f2 2 cos a) < (43) 

(we neglect the term r\jx compared to unit). If rotation is cylin- 
drical and n 2 = then instability occurs if the angular velocity 
decreases with the cylindrical radius. 

Stability properties can be more complicated, however, if 
rotation departs from a simple cylindrical law. Consider, for ex- 
ample, the so-called shellular rotation when is approximately 
constant on spherical shells, SI = Q(r), r is the spherical ra- 
dius. Such rotation can be a reasonable approximation, for ex- 
ample, in the radiative zones of massive main sequence stars 
(see, e.g., Zahn 1992, Urpin, Shalybkov & Spruit 1996). For a 
shellular rotation, we have f2 2 = f2 2 sin 9 and Q 2 = Q 2 cos 9 
where ft 2 = dn 2 /dr. Then, the condition (43) reads 



f2 2 sin a(sin a sin 9 + cos a cos 9) < 0. 



(44) 



If the angular velocity decreases with the spherical radius, 
ft 2 < 0, the instability occurs in the region where 



sinacos(a — 9) > 0. 



(45) 



There is no generally accepted point of view regarding the ge- 
ometry and strength of the magnetic field inside stars. For the 
purpose of illustration, we consider therefore the simplest mag- 
netic configuration that can exist in the radiative zone 



B r = /(r)(l - 3 cos 2 9), B e = F{r) sin0cos6>, 



(46) 



where / and F are functions of the spherical radius alone and 
satisfies the divergence-free condition 



F = 



1 d{r 2 f) 
r dr 



(47) 



The magnetic field (46) has the same geometry as the velocity 
pattern of meridional circulation caused by shellular rotation 
with O = O(r) (Zahn 1992, Urpin, Shalybkov & Spruit 1996, 
Talon et al. 1997)). Therefore, the assumption (46) is a good 



approximation it the held is rrozen-in a slow circulation How. 
Generally, the held in radiative zones can have a more com- 
plex geometry but the simple model (46) allows to understand 
qualitatively how the magnetic held influences the region of 
instability. 

We have for the magnetic held (46) 



sin a = 



cos a = 



e s B 


sin 9 j 


B p 


Bp 


e z B 


cos 9 


Bp 


Bp 



[/(l-3cos 2 <9) + Fcos 2 9], 
[/(l - 3 cos 2 (9) -F sin 2 <9]. (48) 



Substituting these expressions into equation (44), we obtain 



(1-3 cos 2 9) [1 - (1 - q) cos 2 9} > 0, 



(49) 



where q = din //din r. The criterion (49) is fulfilled at 
cos 8 ~ 0, and the region near the equator is secularly unstable. 
Near the rotation axis where cos Sssl, equation (49) yields 



-f <0, 
ar 



(50) 



and this region can be unstable only if the radial held compo- 
nent decreases to the surface. The regions of instability as a 
function of the parameter q are shown in Table 1 . 

There is a principal difference between the cases q < 
and q > 0. The hrst case can correspond to the field generated 
somewhere in the central region of the star, for example, in the 
convective core that exists in massive stars. The second case is 
likely more relevant to the field generated in the surface layers. 
For instance, this can be typical for the late-type stars which 
have outer convective zones where turbulent dynamo operates. 
If the held increases with r and q > then the instability oc- 
curs only in the region with openning angle m 36° around the 
equator. A wide region around the rotation axis with 9 < 54° 
is stable. On the contrary, if the held decreases to the surface 
and q < then there are two different unstable regions. For a 
slowly decreasing held, > q > — 2, the instability occurs in 
a relatively wide unstable region around the equator (with the 
openning angle 9 > 54°) and in a narrow region around the 
axis (with the openning angle arccos 1/y/l — q). If the field 
decreases rapidly with r and q < — 2, then the instability oc- 
curs in a region with 9 < 36° around the rotation axis and in 
a narrow region around the equator (with the openning angle 
arccos 1/^1 — q). Note that this openning angle can be very 
small if the held decreases rapidly. 

4. The growth rate of instability 

Consider the growth rate of the modes assuming that the dissi- 
pative frequences u v , u v and lo x are smaller than the dynami- 
cal frequences u>a, u g and £1 because this case is of particular 
interest for stellar radiative zones. Under this assumption, the 
dispersion equation (15) can be solved by making use of the 
perturbation procedure. We expand 7 as 7 = 7^ + + ... 
where and 7W are terms of the zeroth and hrst order in 
dissipative frequencies, respectively. The corresponding expan- 
sion should be made for the coefficients of equation (15) as 



well. In the zeroth order when dissipation is neglected, equa- 
tion (15) reduces to a quadratic equation, 



7 4 + (2co 2 + u? g + Q 2 ) 7 2 + uo 2 A (u? A + uo 2 g + F ) = 0, 
This equation describes four modes with the frequencies 



(51) 
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1^+40^, (52) 



) 2 ) 2 k 2 
l^+4u; 2 Q 2 |f. (53) 



The hfth root of equation (15) is vanishing in the zeroth ap- 
proximation, 75°' = 0. Equation (51) has unstable solutions 
with Re 7 > if one of the conditions 



2ivj 



f<0, lu 2 a +lu 2 + F 2 <0, 



(54) 



is satisfied. The hrst of these inequalities describes the well- 
known Rayleigh instability, and the second condition corre- 
sponds to the magnetorotational instability (Velikhov 1959, 
Chandrasekhar 1960, Balbus & Hawley 1991). Most likely, nei- 
ther of the inequalities (54) is fulfilled in radiative zones be- 
cause of a strong stabilizing effect of stratihcation. The only 
exception are pertutbations with k || G for which ui 2 = 
(Urpin 1996). Note that, the stabilizing effect of stratification 
is much suppressed for the instability considered in this paper 
(see equation (18)) compared to the Rayleigh and magnetoro- 
tational instabilities, 

The hfth root of equation (15) describes a secular mode 
that is linear in dissipative effects. To calculate this root we 
should keep in equation (15) only the terms linear in dissipative 
frequencies assuming that 75 is linear as well. Then, the hfth 
root is approximately given by 



7s 



(i) 



— Uly 



F 2 +u\ + (77/xV 



F 2 



J A T 

The order of magnitude estimate of this root reads 



(i) 



n 



Wx „ 4 x lO^V^nAg 2 ^ 2 s-\ 



(55) 



(56) 



where X7 = x/10 7 cm 2 /s and V7 = y/10 7 cm/s. Then, the 
growth time of instability is 



r = 1/75 



(i) 



10 2 x7 1 w 2 3 i? 2 1 A 2 yf 2 yrs, 



(57) 



The growth time can be rather short if the star rotates rapidly, 
V 7 ~ 1. 

5. Discussion 

We have considered the diffusive instability of differentially 
rotating magnetized stellar radiative zones. The instability 
considered is a generalization of the well-known Goldreich- 
Schubert-Fricke instability for the case of magnetic stars. In 
non-magnetic stars, this instability occurs if the angular veloc- 
ity depends on the vertical coordinate, O = il(z). The mag- 
netic held, however, changes drastically the stability properties 
of differentially rotating stars. Even in a weak held satisfying 



The dependence of the region of instability on q for a rapid shellular rotation. 
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the inequality (28), the magnetic analogue of the Goldreich- 
Schubert-Fricke instability may occur in a region where £1 de- 
pends on the cylindrical radius alone, £1 = Q(s). The necessary 
condition of instability in this case is a decrease of the angular 
velocity with the cylindrical radius. If dVl/ds < than the in- 
stability may occur in both slowly and rapidly rotating stars. 
If rotation is more complicated, £1 = Q(s, z), then the verti- 
cal dependence of £1 can provide either stabilizing or desta- 
bilizing effect depending on the value and particular shape of 
£l(s, z). In rapidly rotating stars with the equatorial velocity 
V > 10i?iiWo3 km/s, the vertical dependence of £1 plays a 
destabilizing role, and the radiative zones of such stars should 
be unstable. In slowly rotating stars, the vertical dependence 
destabilizes differential rotation if <9£! jdz > and stabilizes if 
dil/dz < 0. In such stars, generally, the instability can occur 
only in some regions near the rotation axis or equator depend- 
ing on the particular shape of £7(s, z). 

In a strong magnetic field with B p > B cr (but still satisfy- 
ing the condition c s > ca), the stability properties are essen- 
tially different. The secular instability does not occur in slowly 
rotating stars with the equatorial velocity V < 10i?nWo3 km/s. 
However, if rotation is rapid and satisfies the condition (42) 
then the radiative zone or its fraction can be subject to instabil- 
ity. The geometry of unstable regions is sensitive to the partic- 
ular shape of £l(s, z). Even in our very simplified model where 
the axisymmetric magnetic field is considered, the instability 
may occur not in the whole radiative zone but only within some 
regions restricted by cones in the meridional plane. For exam- 
ple, in the case of a shellular rotation when il depends on the 
spherical radius alone, the polar angle corresponding to the un- 
stable regions is determined by the radial dependence of the 
magnetic field and is shown in Table 1 . If the radial magnetic 
field increases to the surface then the instability arises only in 
the region around the equator. In the opposite case when B r 
decreases to the surface, the instability occurs both around the 
equator and rotation axis. The openning angle of unstable re- 
gions depends on how fast the magnetic field decreases with 
r. 

The growth time of instabilities is relatively short and, 
likely, it operates in a non-linear regime. We can estimate 
the saturation velocity using the mixing-length model (e.g., 
Schwarzschild 1958) that assumes that the turn-over time of 
turbulence generated by instability is of the order of the growth 
time of this instability. In a vanishing magnetic field, this sim- 
ple approach provides a quantitatively correct estimate of the 
saturation velocity generated by the Goldreich-Schubert-Fricke 
instability in stellar conditions (Korycansky 1991) and accre- 
tion disks (see Urpin 2003c and Arlt & Urpin 2004). Then, the 



saturation velocity in a turbulent cell with the lengthscale A can 
be estimated as t"r(A) ~ A 7(A) or, using equation (57), as 

vt-OMxtu^R^^V? cm/s. (58) 

It is worth remarking that vt is larger for small scale motions. 

Turbulent motions can enhance transport and mix the ma- 
terial in the unstable regions. The coefficient of turbulent diffu- 
sion, vt, can be estimated as a product of the turbulent velocity 
and lengthscale. Then, we have 

VT ~ VT \ ~ 4 x 10 6 X7Wr7 3 2 ^n 2 v 7 cm 2 /s. (59) 

Turbulent diffusion caused by the instability is sufficient to mix 
the unstable zone on a relatively short timescale. Using equa- 
tion (59), we can estimate the timescale of mixing as 

R 2 

W 6 x lOV^MiV 2 yrs (60) 

vt 

that can be substantially shorter than the lifetime of many stars. 
Therefore, elements from a deep interior can likely be trans- 
ported to the stellar surface by turbulence within the unstable 
zones. 

This conclusion is of particular importance for strongly 
magnetized stars with the rotation velocity satisfying the con- 
dition (42). In such stars, the instability leads to a formation 
of non-spherical "turbulent zones" mixing the material of deep 
interiors and the surface layer. In our very idealized model with 
the axisymmetric magnetic field and shellular rotation, the in- 
stability can form two "mixing zones": the cone turbulent re- 
gion around the rotation axis and the turbulent belt around 
the equator. If rotation or the magnetic field are more com- 
plex then the geometry of turbulent "mixing zones" can be 
more complicated. For example, if the magnetic field is non- 
axisymmetric then the distribution of unstable regions should 
be non-axisymmetric as well. The material should be well 
mixed in the "mixing zones" and, likely, this mixed material 
can reach even the atmosphere if the unstable region extends 
to the surface. If the magnetic field is axisymmetric then the 
pecular regions, where the "mixing zones" reach the surface, 
have likely the belt-like structure or form the polar spot with a 
pecular composition. If the field is non-axisymmetric then the 
instability can lead to a formation of the spot-like structure with 
chemical pecularities. 
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